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ABSTRACT 


Fluxes of momentum, heat and water vapour above 
Lake Wabamum, Alberta are dealt with in this investigation. 
Two power-generating plants discharge heated water from 
condenser-cooling systems into the lake and the fluxes 
above this heated water are the main concern. Measurements 
of the fluxes by fast-response instruments are presented 
for two warm days in August and one cold day in October. 
Temperature, dew-point and wind profiles are included for 
the October date only. 

The August fluxes are calculated from one day's 
data collected above the heated water and from a similar 
day"s data collected above the undisturbed lake. It is 
shown that the fluxes of heat and water vapour above the 
heated water were approximately 5 times greater than those 
above the natural lake. 

Et is noted. that free convection conditions exist 
above the heated water for the October situation because of 
a 15°C difference between air and water temperature. The 
October profile and fast-response instrument estimates 
of fluxes are compared and found to agree well for sensible 
heat flux, less so for vapour flux and not at all for 
momentum flux. An upward momentum flux as measured by the 


fast-response shear-stress meter is discussed and is shown 
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to be linked with a large-scale turbulence upon which 
the normal small-scale turbulence is superimposed. The 
possibility of the large-scale turbulence being associated 
with organized convective eddies is pointed out. 

A liquid water flux associated with a dense fog 
during the October data collection is considered and 
found to be an order of magnitude less than the water vapour 
na Hy bees 

Sensible and latent heat fluxes due to an observed 
positive mean vertical velocity in October are calculated 
and discussed. Each is shown to be significant but less 


than the corresponding turbulent flux. 
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CHAPTER I 
INTRODUCTION 


Over ine past few years in Central Alberta, 
considerable scientific and community interest has been 
focussed on a body of water about forty miles west of 
Edmonton. This lake, Lake Wabamun, is both large enough, 
(three by twelve miles) and close enough (42 miles) to a 
large urban population centre to serve as a well-used 
recreation area. Thus, summer cottages, beaches, and 
other similar facilities including a provincial park border 
much of the lake. In addition, accessible coal has been 
found on the shores of the lake making the area an ideal 
location for thermal power plants. Coal for fuel and water 
for steam and condensing purposes are convenient and 
abundant. As a result, Calgary Power has established two 
génerators’)(599 and: 300 MW) on‘oppositershores of the’ lake 
and these generators discharge water from their cooling 
systems back into the lake. It has been claimed that the 
resultant plumes of heated water have caused increased weed 
growth and other biological and zoological changes 
detrimental to the recreational aspects of the lake. 
Consequently biologists, zoologists and Givil engineers 


from the University of Alberta became involved in determining 
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the validity of such claims through the Wabamun Lake 
Thermal Project (e.g. see Nursall and Gallup, 1971). During 
1972 the Meteorology Division of the University entered the 
investigation in an attempt to examine the fluxes of 
momentum, heat and water substance from both the heated 
plumes and the ambient lake surface. I was fortunate 
enough to take part in this micro-meteorological project 
and the following chapters deal with the techniques of this 
project, and some of the results. 

Although the bulk of boundary-layer flux measurements 
have been made over land, the problem of estimating heat 
and water exchanges between air and a water surface has 
received woridwide attention for many decades. At first 
most estimates of evaporation were obtained from empirical 
flux-gradient relationships involving wind speed, surface 
water temperature and air temperature and moisture content. 
These formulae were frequently used in water-level studies 
of lakes and reservoirs and similar hydrologic im tigations. 

A logical extension of this flux-gradient : ationship 
assumption leads to estimates of all the boundary-layer fluxes. 
The technique is simply to apply a mixing-length type approach, 
where the flux of a quantity is given by a constant times 
the gradient of the concentration of the quantity. <jAsy is 
pointed out by Robinson (1966), this approach could be 


extremely useful as it relates phenomena physically compounded 
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of micro-scale interactions with standard meteorological 
observations. Thus, using synoptic-scale data only, 
estimates of the Bee laa: fluxes are available, 

e.g. for input into dynamic models. Such a method is very 
approximate, however, and the resultant flux estimates 
approach the real fluxes only when averages over vast areas 
of ocean are considered. The same method can be applied 

to lakes that exhibit horizontal homogeneity in atmospheric 
moisture and temperature and in surface water temperature. 
Edge-effects complicate the procedure, but for sufficiently 
long time intervals, fairly accurate estimates of the fluxes 
can be obtained if the above restrictions are met. 

One of the most comprehensive studies of evaporation 
from a lake was the well-known Lake Hefner project (e.g. U.S. Dept. 
PDLALLOy, 1954). The purpose of the investigation was to 
obtain formulae to predict evaporation from Lake Mead before 
Lake Mead existed. Because inflow and outflow for Lake 
Hefner were well known, a water-budget method was used as a 
control estimate of evaporation. Observations from a meso- 
scale network around the lake were used to compute evaporation 
estimates using an energy-budget technique as well as various 
empirical flux-gradient procedures similar to those mentioned 
previously. Generally successful comparisons were then made 
between these estimates, and the water-budget estimates, of 


average monthly evaporation. 
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As the Monin-Obukhov Similarity Theory (e.g. Lumley 

and Panofsky; 1964) gained rapid, widespread acceptance, | 
and with improved instrumentation and computing facilities 
profile methods became more sophisticated, while direct 
measurement of the turbulent fluxes became possible. 
Similarity theory provided a sound theoretical basis for 
the form of the profile-flux relations and also introduced 
the: inclusion of Stability effects. Thus, estimates of flux 
averaged over minutes, instead of days, were feasible. 
Additionally, the last fifteen years has seen increased 
use of instruments with responses fast enough to measure all 
Significant turbulence fluctuations. If the responses of the 
instruments are properly matched, then the covariance of two 
measurements is directly proportional to the appropriate flux. 
peti axonks temperature and wind sensors have been used 
for some time to measure heat and momentum flux. More 
recently, the use of the Lyman-alpha humidiometer has 
permitted direct measurement of turbulent water-vapour flux. 
The disadvantages of these methods are the cost of the 
instruments, and the large Bevin of data and data anlaysis 
required for even short-period calculations. The data 
collection and analysis problems become practicable only 
with the use of modern computers. These are hardly the 
methods for calculating evaporation on a long-term or 


large-scale basis, but they are the only methods which give 
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information about the structure of the boundary layer 
itself. 

Two of the largest expeditions to use the improved 
profile methods, and direct-measurement approach to the 
problem of evaluating fluxes over water, were the Indian 
Ocean Expedition (Badgley, Paulson, and Miyake, 1968) and the 
Barbados Oceanogrophic.and Meteorological Expedition (BOMEX) 
(e.g. Holland 1972). Both projects were concerned with 
determining fluxes in the atmosphere above large expanses 
of fairly warm ocean. BOMEX, in particular, attempted to 
evaltiate' the fluxes using every pdéssible “technique for 
comparison purposes. Thus, in addition to the two previously- 
mentioned methods, budget methcds, and a procedure based on 
Kolmogoroff's Law were also used (Pond et al. 1971). Budget 
methods siiptor conservation equations for heat, water vapour, 
or horizontal velocity to express boundary-layer fluxes as 
a function of observable quantities. The fourth procedure 
involves calculation of the turbulent-transfer coefficients 
and vertical fluxes through the use of Kolmogoroff's Law, 
which relates the energy eoedtien and wave number in the 
inertial subrange. The observations used for the direct 
covariance calculation of fluxes can also be used in this 
method. 

Both profile and fast-response type measurements 


were made recently at Lough Neagh, a lake in Northern Ireland, 
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by Sheppard, et al. (1972). In fact, the Lake Wabamun 
experiment was similar in design to the Lough Neagh 
operation, though far less extensive in scope than the latter. 
The combination of profile and fast-response measurements, 
allowed both fhohss to make use of three of the four 
independent methods of determining fluxes in the atmospheric 
boundary layer. 

As was mentioned previously, the empirical flux- 
gradient relationships described earlier are intended for 
use over large bodies of water that exhibit horizontal 
homogeneity. This is hardly the situation at Lake Wabamun, 
where it is the departures from horizontal homogeneity which 
we wish to examine. The water-budget method, used on Lake 
Mead, provides only an estimate of evaporation averaged 
over the entire lake for fairly long time periods. For 
this reason, and because of the probable existence of a large 
underground source for Lake Wabamun, it was deemed to he 
unsuitable. . Thus, the only techniques appropriate to the 
Wabamun project were those supplying information about the 
vertical fluxes at some point Aus the atmosphere. Of the 
four independent methods left, three can be utilized if 
profile and fast-response measurements are available. With 
this in mind, these were the two types of measurements made 
at Lake Wabamun during 1972. 

Useful data were obtained on three days, two of which 


were warm days in August, while the third, an-october, had 
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7 
air temperatures near freezing. Different wind directions 
permitted measurement of fluxes above both the heated 7 
plume and the undisturbed lake surface for the August 
Situation, but only the fluxes above the plume were obtained 
in October. It is the intent of this thesis to compare 
various methods of estimating fluxes at some point in the 
atmosphere, as well as to compare fluxes from the heated and 


unheated water of the lake. 
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CHAPTER II 
COLLECTION OF THE DATA 
2.1 The Observing Site 


Lake Wabamun is situated in gently rolling, 
forested terrain, 67 km west of Edmonton, Alberta. The 
lake lies in an east-west direction, being approximately 
19 km long and on the average 4.3 km wide. A more 
complete list of the physical characteristics can be 
found in Table 1. The Wabamun power plant, the larger 
of the two generating stations, is on the north shore of 
the lake near the eastern end. The eastern tip of Point 
Alison was selected as the observing site for the 1972 
experiments in order to take advantage of air flow over 
the heated water discharged from this station (see 
Figure 1). The point protrudes into the lake between the 
power plant's intake canal to the west, and the discharge 
canal to the east. The only acceptable wind directions 
at this site were east to southeast with an over-water fetch 
pEnee less than 2.2 km, and west to southwest with an 
overwater fetch of not less than 5.6 km. East to southeast 
winds carry warm water from the outlet westward around 
Point Alison and the air reaching the observing site has 


a fetch of at least 0.5 km over the heated water plume. 
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TABLE 1 


Morphometry of Lake Wabamun 


Elevation 

Area 

Volume 

Length 

Maximum breadth 
Mean breadth 
Maximum depth 
Mean depth 
Shoreline length 


Area of Surface drainage 
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Figure 1. - Lake Wabamun, Alberta - showing experimental site at 
Point Alison. 
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Figure 2. Point Alison and vicinity, Lake Wabamun, Alberta. 
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1 
In a southwest wind the warm water plume is carried east 
of the point and air reaching the site has a trajectory 
over the normal lake surface only. 

An enlarged view of Point Alison showing the 
experimental site is shown in Figure 2. A trailer, housing 
recording and data processing equipment, was located near 
the water iine, and an 8 m tower was erected in water of 
cepth 0.59m abouty25 a out Erom: the lakeshore. “Point 
Alison itself is heavily forested and rises to a maximum 
elevation of 18 m approximately 0.6 km westnorthwest of the 
tower site. As is shown by Figure 25» Che £Fip Gs, Hie. point 


may cause flow interference in a wind from due east. 
2.2 Instrumentation 
2s 2k Slow-Response Instruments 


The slow-response instruments listed in Table 2 
were used to obtain average values of wind, temperature, 
humidity, and net radiation for time periods of fifteen 
minutes or more. The wind speed, air temperature and air 
dew-point instruments were mounted at the four listed 
heights on An 8 m.tower *to provide*the profile data. The 
non-logarithmic heights were necessitated by the use of 
a tower with fixed booms, designed for use on land, but 
sitting in 0.5 m of water at Lake Wabamun. Net radiation 


and wind direction instruments were located on separate 
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masts away from the tower and their signals were recorded 
in analog form on strip charts. Dew-point temperatures 
and water temperatures were read and recorded periodically 
from meters, while Assman psychrometer readings were taken 
at irregular intervals at low heights on the tower as 
additional independent checks of air temperature and 
moisture. 

The temperature profile utilized linear thermistors 
and a set of operational amplifiers. The thermistors were 
mounted in reentrant shielded tubes which were ventilated. 
The output was arranged to give the ambient temperature at 
One level, and differences from that temperature at the 
other levels. Maximum sensitivity to the temperature 
gradient can be obtained in this way. Data readouts were 
accomplished by driving integrating motors with number 
wheels, which were photographed every fifteen minutes to 
give a mean temperature. 

The dew-point profile instrument had four sensing 
heads. They were connected to the thermistor elements of 
the temperature profile, and ventilated together by a small 
suction pump. The dew-point detection was accomplished by 
sensing light scattered from condensate on a mirror which 
was cooled to dew-point by a thermoelectric cooler (Peltier 
effect). To obtain profiles, the four elements had to be 


switched on manually in turn, and read in sequence. 
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14 
The sensitive cup anemometers mounted on the end 
of each of the tower booms measured total horizontal wind 
at each level. The cup generators operated counters 
located in the trailer which were photographed every fifteen 
minutes, permitting the calculation of mean wind speeds 


for periods of a quarter hour or more. 


Bene Fast-Response Instruments 


Fast-response instruments at the tower site are 
listed in Table 3. All channels were sampled’ at 100 
readings per second and filtered at. a (Cut-off Eregquency ‘of 
20 Hz before conversion from analog to digital forn. 
Analog-to-digital conversion, gross error checks, and 
recording of; up| to 8 channels of data on digital magnetic 
tape were carried out in the field following procedures 
described by Honsaker, McDougall, and Oracheski (1972). 

A shear-stress meter, consisting of three Gill 
propeller anemometers mounted on orthogonal axes, was used 
for the measurement of eddy flux of momentum. It was 
mounted on a separate mast so that it could be moved 
independently of the main tower,thus avoiding wind inter- 
‘ference effects. 

The sonic anomemeter-thermometer was a continuous- 
wave model similar to that described by Kaimal and Businger 


(1963). The instrument was capable of measuring fluctuations 
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16 
in the vertical wind component and temperature simultaneously. 
It was mounted on the main tower 3 m above the surface aoe 
water. 

A Lyman-alpha Humidiometer was attached to the 
same mounting as the sonic anemometer-thermometer. The 
Humidiometer produces a signal from the absorption of 
ultraviolet light around the Lyman~alpha wavelength 
(121.5 nm). In air, Lyman-alpha absorption is almost 
entirely due to water vapour, and thus provides a measure of 
absolute humidity. The strength of the absorption permits 
the measurement to be made over a short sampling path 
(15 mm) which in turn means that the instrument is capable 
of very fast response. It was used in conjunction with the 
sonic anemometer to measure the eddy flux of latent heat. 
Uerieerare resistance thermometer was mounted 
on its own mast with a Gill propeller anemometer measuring 
vertical wind either on the same mast or on the adjacent 
shear-stress meter mast. The resulting simultaneous 
measurements of temperature and vertical wind provided an 
estimate of the sensible Heat cant which was independent of 


the sonic anemometer-thermometer estimate. 
2.2.3 ° Instrumental Array 


The arrangement of instrumentation at the site is 


shown in Figure 3, a photograph of the tower looking southward 
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18 
across the lake. The profile instruments were located near 
the ends of 1m booms oriented southeastward toward the 
hot water plume. Thus the booms were aligned approximately 
into the wind during data collection and blockage effects 
were minimized, The sonic anemometer-thermometer and the 
Lyman-alpha humidiometer were attached to a boom extending 
southwestward from the tower 3 m above the lake surface. 

The shear-stress meter, resistance thermometer, net radiometer, 
and wind direction indicator were movable and were arrayed 
across wind from the tower before each day's runs. All 

masts and instruments were aligned with the vertical by using 


transit readings at three points along the shore. 
nwo “Calibrations 


Considerable effort was expended calibrating the 
instruments. The anemometers (cup and propeller) were 
calibrated against a hot-wire anemometer and a pitot tube 
in the University of Alberta Civil Engineering Department's 
wind tunnel, and against a secondary standard anemometer in 
a wind tunnel of the Defence Research Establishment at 
Suffield, Alberta. Further calibrations of the cup anemometers 
in both the laboratory and the field were necessary for low 


wind speeds (2 m/sec or less) because of different responses 


among the anemometers in this range. 
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The thermistor elements were calibrated using a 
temperature-controlled air flow in the laboratory spanning 
the entire range of temperatures encountered in the field 
experiments. Temperature dependent noise in the signal, 
discovered during these calibrations and affecting the 
integrating motors, eventually necessitated the discarding 
of the temperature profile data from the two August data 
collection days. The dew-point hygrometer was calibrated 
in both field and laboratory environments. As well the 
Assman psychrometer was used as a check on air temperature 
and moisture content during the actual data collection. 

The fine-wire resistance thermometer was calibrated 
simply by measuring its resistance at several temperatures 
and’ fitting a, linear relation to, the points... The 
humidiometer calibration was determined by taking the reading 
for dry air. Because of the exponential relation between 
absorption and the water vapour amount, one reading defined 
a unique curve and complete calibration. 

The sonic anemometer-thermometer has a built-in, 
electronically simulated, 30° phase shift. This, phase, shitt 
was recorded both manually and directly onto magnetic tape 
permitting calibration to be carried out in the field before 


or after data collection 
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2.4 Water Temperature Surveys 


Temperature measurements in the lake were carried 
out from a small boat using a calibrated thermistor read 
through a YSI Model 42SL telethermometer (See Table 2). 
Surface temperatures were recorded by holding the thermistor 
in the upper layer of water with about 2 inches submersion. 
Temperature profiles were obtained by anchoring the boat 
and then lowering the thermistor by 2 foot increments. 
Bottom temperatures were recorded by dropping the weighted 
thermistor head into the mud bottom and waiting unti. 
equilibrium was obtained. 

In August the boat's position was found by means of 
two simultaneous theodolite bearings from instruments set 
up Om shore about 1200 feet apart. The base line distance 
was measured to within 0.1 ft and the angles recorded to the 
nearest 5 seconds of arc. Citizen's band radios were used 
for communication between the boat and theodolites,. 

Dense steam fog from the lake made it impossible to 
use the same procedure to determine the boat. Ss position in 
October. Because of the extremely low visibility, only one 
traverse of the plume directly upwind of the tower was 
possible. The procedure was simply to start at the tower 
on a heading directly into the wind, continuing into the 
lake on the same heading until homogeneous lake temperatures 


were obtained. The boat's speed was estimated by measuring 
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24 
the time for objects on the water to pass from bow to stern 
of the boat in a manner similar to the casting of the log 
on a Sailing ship. Of necessity then, the distribution of 
upwind water temperatures was subject to considerable 
uncertainty. However, because the width of the plume matched 
the known width on days with similar winds, the data for 


October were accepted. 
2.5... LOCAL Weather Conditions 


BOed AUCUSst 24 and 25, L972 were warm (7/5-"Fy 24° C) 
summer days with only eS tone cloud, “OnvAuguse 24 - Tight 
(3-4 m/sec), steady southeast winds prevailed through the 
day and into the evening. As was discovered during previous 
water temperature surveys with similar winds, the hot water 
plume was ee westward around Point Alison and directly 
on shore at the tower site (see Figure 4). Surface water 
temperatures reached 30° C in the plume so that the air 
arriving at the tower had been receiving heat from the 
lake surface for some 600 m and hence was likely to be 
unstable. The plume remited along the Point Alison shore 
until after data collection had ended. Data were collected 
from 1330 to 2045 MDT (Mountain Daylight Time). 

A weak trough passed through Alberta during the 
night of August 24 and winds were steady from the southwest 


at about 2 m/sec through much of August 25. These winds 
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22 
carried the very shallow layer of hot water east of Point 
Alison and air reaching the tower site had passed over 
the cooler (22°C), ambient lake surface (See Figure 5). 

This meant that neutral or slightly stable conditions in 

the surface boundary layer could be expected to prevail 
through much of the day's data collection period. In the 
evening, wind speeds dropped to near calm and wind direction 
became variable. These two days were ideal for comparison 
of fluxes over the plume with fluxes over the rest of the 
lake. Useful data collection occurred between 1200 and 

1630 MDT on August 25. 

On the morning of October 21, 1972, a low pressure 
area in the southern parts of British Columbia and Alberta 
produced overcast stratocumulus cloud conditions, some periods 
of snow, and east winds in the Lake ore region. By 
mid-afternoon the winds had veered to eastsoutheast and this 
direction was maintained well into the night. Once again 
these winds brought the heated water in along the shore of 
Point Alison in a situation similar to that shown for 
August 24 in Figure 4. Thus air with a temperature near 0° C 
had a trajectory of at least 300 m over water whose surface 
temperature approached 17° C (See Figure 6). This temperature 
gradient produced unstable conditions in the atmospheric 
surface boundary layer. A steam fog formed over the entire 


lake and was particularly dense over the warm waters of the 
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CHAPTER III 
HBO eCAT ANA IDBIGEES 
Ser - Introduction 


In the past twenty-five years investigators of 
the surface boundary layer have used three main approaches 
to study and describe the motions of the air near the 
ground. The earliest method examined the average profiles 
of wind and temperature and through existing theory tried 
to infer turbulent structures and flux magnitudes. For 
the sake of simplicity the bulk of this theory was developed 
for an infinite, homogeneous, level surface. A large lake 
approximates these conditions very well and the profile 
approach was chosen as an appropriate one to apply to the 
Lake Wabamun data. The most serious departure from the 
profile-theory assumption at Lake Wabamun was the change 
in surface temperature because of the heated water but it 
was believed that the profiles would have adjusted 
sufficiently at the levels of interest for this not to be 
a serious problem. 

Another approach usesaccurate, fast-response 
instruments to examine the structure of the turbulence and 
the fluxes etc. directly. This method was also employed at 


Lake Wabamun. The third approach in common use involves 
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25 
numerical modelling of the surface boundary layer. The 
computer has proved particularly useful in extending theory 
to cover more complicated flows such as occur over changes 
in surface roughness and temperature, problems that have 
received increasing attention over the past fifteen years. 
The results of one such model will be discussed very briefly 
in chapter four. The theory of the other two approaches 
will be developed in this chapter. 

The *imtent lof *this chapter is to provide a logical 
development of the theoretical formulas used with the Lake 
Wabamun data. A more complete treatment of surface boundary- 
layer profiles and turbulence can be found in a number of 
texts (see e.g. Lumley and Panofsky, 1964) and their contents 


need not be duplicated here. 
3.2 Profile Theories 
Jeeerils PLLNGEOdUc tLon 


Most of the recent advances in understanding of the 
distribution of wind, temperature and moisture near the 
earth's surface have been brought about by the use of the 
Monin-Obukhov "similarity theory." Because of its importance, 
and because much of the present surface boundary-layer 
profile theory is based on it, similarity theory has been 


chosen as the starting point for this discussion. 
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The complexity of the problem makes it impractical 
to completely and theoretically edn the state of the 
surface boundary-layer for the perfectly general case of air 
moving over a varying surface. For this reason numerous 
assumptions were found to be necessary in order to simplify 
the theoretical description of the mean flow. 

When applied to meteorology, the equation for the 
mean motion in vector notation can be written as (see e.g. 
Lumley and Panofsky, 1964): 
Ba aN Be ae SARE AT. 8M, ie D 
d 7? ? 


=) 
Here V is the horizontal wind vector, f£ the Coriolis 


_> 


ecu taaell the density, Vo the pressure gradient, Y the 


5 


ch 
Vv 


tangential stress vector acting on a horizontal surface, 
and k Sea aakerel unit vector. For micrometeorological 
studies of the surface boundary layer the magnitudes of the 
Coriolis force and the pressure gradient force are small 
enough to be neglected leaving as the expression for the 
mean flow: 
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Two standard surface boundary-layer assumptions are that there 
is no external pressure gradient and the earth's rotation 
need not be considered. Both assumptions limit the scales to 


which the results of theory are applicable. A condition 
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27 
necessary for the theoretical development is that the 
turbulence characteristics exhibit stationarity ina 
statistical sense which implies that the mean flow is in 
equilibrium. Thus stationarity becomes a third assumption. 
A further restriction which limits the applicability of the 
results, reduces the problem to a two dimensional one by 
assuming horizontal homogeneity and aligning the horizontal 
axis with the mean flow. Finally, molecular transport and 
radiative flux divergence are neglected although the latter 
assumption may not be realistic in stable conditions. 

Under these restrictions the shear stress or 
momentum flux and the heat flux vary little with height near 
the graund. In fact the surface boundary layer is frequently 
called the "constant-flux" layer defined as the layer in 
which the shear stress varies by less than 20%. Stress 
divided by density has units of velocity squared so that the 
constant shear stress condition makes it possible to define 


* * 
a "£rickion velocity yay, constant with height ,<by: 


% 
wos fe 
Pp (31) 


where Y is the surface shearing stress oud is density. 
For the turbulent flow problem simplified by the 
oreviously-mentioned assumptions, Monin and Obukhov (1954) 
argued that the terms governing the flow characteristics, 
andhence the important terms for dimensional analysis 
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acceleration of gravity, T a mean air temperature, Va a 
mean air density, H the heat flux, and z the height above 
ground. This set of five variables involves four dimensions, 
mass, length, time and temperature, which, by the 
Buckingham -‘7{ theorem, means that all mean flow parameters 
are expressible as a function of one nondimensional quantity. 


This quantity was chosen to be z/L where 


ees a Lab Lp ee Derry Ki 
h G/T) (H/ga) ka H ee 


The constant c,is the specific heat of air at constant 
pressures’ Thesvons Karman constant k is included arbitrarily 
because it is dimensionless. This constant was used 
historically to describe the wind profile and has a 
numerical value of about 0.4. All mean turbulence character- 
istics, when placed in nondimensional form, should be functions 
OL i2/iss 

A scaling length L and a scaling velocity a have been 
found for the flow and for completeness a scaling temperature 


* Fe . . . 
T has been determined from the same dimensionally-important 


quantities (see e.g. Lumley and Panofsky, 1964): 
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The wind shear is related to the stress by 


= - 
pat be ser also) Y (3-4) 
Va 02 

wWhere*KE 25 the scalar eddy Viscosity: *Since*K .*is"a 


M 


scalar, the stress and wind shear have the same direction 
everywhere in the surface boundary layer. But the stress 
is constant in this layer so that the wind shear is constant 
and the wind direction remains constant with height and 
equals the stress direction. Because only one direction is 
involved the vector notation can be dropped and the quantities 
considered as scalars. 

The heat flux H is related to the potential~-temperature 


gradient in a similar manner: 


Hie Saye Ky, 20 (3-5) 


where Ky is the eddy conductivity and @ is the potential 


temperature, 
Ba 2etuerr Oni lLesuanANeutral and Near~Neutral Conditions 


Tt has long been known that.the velocity profile 
in neutral air is logarithmic in height. A simple derivation 


notes that the ratio 
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30 
has'dimensions of length. As well, in neutral air the 
lapse rate of potential temperature is zero and by (3-5) 
the theat flux Wimust be zéerovalso. This means that L 
goes to infinity for neutral stability leaving the height 
z as the only remaining characteristic length. Traditionally 
the von Karman constant is included in the equation so that 
(3-6) is set equal to kz and the differential equation for 


the wind profile becomes 


OV 2 pire: (ea7) 
Oz ; 


Its integrated form is 


Vai is bn taal (3-8) 


The term sz. 1sathe “noughness, length >. and ‘arises from the 
integration as the height at which the wind velocity becomes 
Zero. 

Perfectly neutral conditions in the atmosphere are 
bare, of course: so: that, 2t.1s desirable to extend the theory 
to cover non-neutral conditions by including the effect of 
the buoyant energy production term. 

The term z/L is an indicator of stability based 
upon the heat flux H contained in the denominator Of IE. 

For neutral conditions H = 0 and (z/L) = 0 while for 


potential temperature decreasing with height (i.e. unstable 
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conditions) the heat flux is upward or positive and 2/L 
becomes negative. In stable conditions, the potential 
temperature increases with height, the heat flux is down- 
ward or negative and z/L increases in the positive sense 
with increasing Stabviaty. If z/ligis smarl ain magnitude 
and negative then the rate of production of convective 
energy agi p 4) is positive but small compared to the 
rate of production of mechanical energy (ux? JV/ Oz). 
Under these conditions the boundary layer is said to be ina 
state called forced convection. According to Lumley and 
Panofsky (1964) this occurs when the magnitude of z/L is 
less than 0.03. Small deviations from neutral stability, 
namely forced convection and slightly stable stratification, 
will be considered next in this discussion. 

Following the lead of Monin and Obukhov, a non-: 


dimensional wind shear S can now be introduced as 


Oe ee sR ee oly. (3=9} 


Considering (3-7) it can be seen that this is unity in 
neutral air. A nondimensional temperature lapse rate R can 


be analogously defined by 


RO wee oe Oe (3-10) 
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But similarity theory predicts that all nondimensional 
mean turbulence characteristics should be functions of z/L 


so that (3-9) and (3-10) become 


So i ee eee aye EZ 
u™ dz 
and 
sh afi bien adie an dete Ye (3-12) 
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Thus the understanding of the behavior of QD, (2/b) and 
D(z/t) become essential to the theoretical development 
of expressions for non-neutral profiles. 
As was done by Monin and Obukhov, D, can be 
expanded in a Taylor series. Since the nondimensional 
profile equals unity in neutral conditions, @, (0) equals 


one and for small 2/L 
D, ‘ea die ee A_#. (321-3) 
a 


where re) Ae* ey constant. -integration of *(S-12P)"atter 


substituting (3-13) yields the log + Jinear velocity profile 


Vee Eo + Az| (3-14) 


provided’ 2° >7"2,. Lie constant must be determined from 


observations in near-neutral conditions. 
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The evaluation of the term z/L often poses a problem 
in the experimental situation because L depends upon the 
venticaltheatdflux. If only the profiles are available 2/L 
cannot. be determined. Panofsky (1963) defined a new scale 


length L' that depended Only upon the gradients as 


a 6 ic (JV/Jz2) T (3-15) 
Rg (XO/dz ) 


where L' and L are related by 
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Pa aya We | | (3-16) 


The L' expression is derived simply by substituting in 
(3-2), for fhe and: Hi from, (3-4) vand 13-5). 
The relative size of Ky and Kus has long been a 
topic for debate but almost all agree that the ratio 
K/Ky is close to unity for neutral air, increases in unstable 
air and decreases in stable air. Again following Panofsky 


(1963) it can be postulated that KL /Ky depends only on 


z/L implying that S = G,(27L™. Then “20m smali.2/i" 
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most investigators proposing values between 4 and 6 for the 
range where the nonlinear terms of the Taylor expansion 
can be dropped. 


Now from (3-4) ~ 


a 
Ky ew . Reha (3-18) 
(9VA2) Z, (2/1) 
and from (3-5) 
eee es shld 2 eed watt 6 (3-19) 


b.(2/.) = O (z/.) ky /K, es 


If the ratio K/Ky is near unity in near-neutral 
air, then the temperature profile for forced convection and 


slightly stable conditions is given by 
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Zo 
where ©, is thexpotential, bemperatore’ at a= 25 
The term z/L is not’ the only term used as a measure 
of stability in the surface boundary layer. A logical 
stability parameter is the ratio of the rates of buoyant 
energy production and mechanical energy production. This 
ratio is commonly called the flux Richardson number R. and 


likez/L, is dependent upon the fluxes so that it is 
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35 
advantageous to evaluate it in terms of the gradients. 
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The term Ri has been defined as the gradient Richardson 
number but because it is used so commonly it is usually 
referred to as simply the Richardson number. It is related 


to the other stability parameter z/L as well: 
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The Richardson number has the great advantage that it can 
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be estimated from measurements of winds and temperature 
at two levels. Because of this ease of calculation, it has 


become the standard stability parameter. 
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Sia. sar Profivlescinerree Convection 


Next to be considered is the extremely unstable 
case, the limit of which is called free convection. 
According to Lumley and Panofsky (1964) "free convection is 
the state of the boundary layer in which the vertical transfer 
of heat and momentum by mechanical turbulence can be 
neglected compared to that by heat convection." As Lumley 
and Panofsky have pointed out, this does not mean that the 
quantity of mechanical energy produced is negligible 
compared to the quantity of convective energy. It has 
been determined empirically that free-convection conditions 
exist when the absolute value of the Richardson number is 
Wnty. ‘According, to ($=25) "and: (3-22), this’ means, that 
mechanical energy production and heat energy production are 
of the neil order of magnitude. In fact Pandolfo (1966) and 
Priestley before him have claimed that free convection exists 
for Ri= -0.03. Apparently then mechanical turbulence is 
a very inefficient transporting agent compared to heat 
convection. Following Lumley and Panofsky (1964) again, 
the explanation of this lies in the fact that mechanical 
eddies have small wavelength and are more nearly isotropic 
while convectional eddies are larger and less isotropic, 
providing a much larger correlation between horizontal and 
vertical wind. Thus the small, isotropic mechanical eddies 


transfer relatively little momentum and heat compared with 
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the large convectional eddies which have a directional 
preference in the vereaes, 

The absolute value of all the stability parameters 
is directly dependent upon the height z so that in every 
unstable constant-flux layer there is a height great enough 
CO perme fréee—convection conditions to exist, and a height 
small enough to allow forced convection assumptions to hold. 
The temperature gradient determines the absolute heights 
at which the two extremes of turbulent transport in the 
unstable surface boundary layer can be considered valid. 

Because mechanical transport can be neglected in free 
convection, u* is no longer important in the theoretical 
development. The dimensionally~-important quantities for 
free convection become Se H, and g/6, with dimensions of 
length (L),.mass (M), time (T), and temperature (8). This 
means that, by the Buckingham ~-7f theorem, nondimensional 
arrangements of the variables are equal except for a 
multiplicative constant. Thus the units of the relevant 
variables in the free-convection potential temperature 


profile (JO/dx ) are 
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can be determined from 
DB mes Ma eh) 
dz x ee 5 
rete em lene ih bor 0 uf ienueds 


Solving for the powers to keep the equation dimensionally 


consistent yields « = -4/3, & =e Cx -1/3, i =e 


and the temperature equation becomes 
! J ~4 
5 


SES ez eae cme (= Zz (3-26) 
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where C is a constant. Integration yields 
(é ) 
ee 


Here B is an integration constant. The value of C can 
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be found if heat-flux measurements are available. From 
measurements by Priestley (1959) its value is 1.07 while 
Dyer (1965) has determined a value of C = 0.83. 

The eddy conductivity, defined by aS ay ( yiey oe) 


can be expressed using (3-26) as 
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a - 
Ve Vabad ute ( qH } are (3-29) 
Oz Kg Cp 8, 


There are two problems associated with this result. In 

the original development of free-convection relationships, | 
u* was assumed to be unimportant in determining the profiles 
Dit wi tolappears explicitly in this wind profile, This 
Seriously reduces the credibility of (3-29). In addition 
the variation of the ratio K/Ky with height in free 
convection is unknown. If this. ratio is independent of 
height then the wind profile follows a sae power law 


and becomes after integrating (3-29) 


xo “4 -4 =4 
V (z,)- Ven Ss uU (4 | Zz, rs 2, | (3-30) 
Kay te) 6, 


According to Panofsky, Blackadar and McVehil (1960) for 
example, this fits wind observations on convective days 
with moderate wind speeds quite well. 

Pandolfo (1966) has derived an alternative 
expression for the free-convection wind profile which still 
depends on u* explicitly but no longer makes the dubious 
assumption of a constant Kj /Ky ratio. Pandolfo based his 
derivation on a rearranged form of the free-convection 


temperature-gradient equation (3-26): 
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where g*: “HAdke, p) and C, = 3k *% C. Pandolfo then 


postulated that the free-convection wind profile took the 


similar form 


Ku (3-31) 
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He then presented strong empirical evidence for z/L = Ri 


which, using. (3-25), leads to 
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Then integrating (3-31) after substituting (3-33) allowed 
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Pandolfo to express the free-convection wind profile as a 
Z power law 


ne 
x oe % 3 

Viz,)-Vizy = ue 6 (a) [La] 2 2 | (3-34) 
h 3 i LB 


6 
Pandolfo (1966) then presented evidence that wind profiles 


in free-convection conditions fit (3-34) well. Later in 
this discussion (3-31) and (3-34) will be tested against 


free-convection wind profiles obtained at Lake Wabamun. 
3.2.4 General Profiles in Unstable Air 


Next to be considered are the profiles for unstable 
conditions in between the two extremes of free and forced 
convection. Priestley (1955) has suggested that the 
transition from forced to free convection is quite sharp 
taking place ata) Richardson, number of ~0.03., Other 
measurements, including some by Priestley (1959), indicate 
a gradual transition. Apartetrom this debate; there 1s 
obviously some advantage in being able to describe the 
structure of wind and temperature through the entire unstable 
range by means of one relationship. Such an equation was 
first found by interpolation between free and forced convection 
by Ellison (1957). The same equation was independently obtained 
by Kazanski, Yamamoto, Panofsky and Sellers and the relation- 


ship. has become known as the Keyps equation after the first 
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letter of the names of ‘the originators. This relation 


between the nondimensional wind shear S and z/L is 


eo ey Germs: patie (3-35) 


Lumley and Panofsky (1964) showed that the relation is 
quite good with the constant y =14. 

An analogous equation using L' instead of L, and 
hence using temperature profiles instead of heat flux, is 


obtained by substituting (3-16) .ine(3-35). 


Shiels y (2A!) Satire rd (3-36) 


Here y is defined as ¥ Ki /Kys Lumley and Panofsky (1964) 
claimed that y is also constant at about 18 which implies 
that K/ Ky must also be constant at 1.3. This resuit is 
highly controversial and disagrees with the observation of 
Pandolfo (1966) that K,/K, = 1/ GZ, (Z/L). 
Following Panofsky (1963), the condition 
= Ds (z/) = RE fase) 
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The roughness Z, is always small compared with L' and because 
the integral is well behaved at the origin (3-37) can be 


written 
OD se bit Se laee A by Vii hag I (3-38) 
h Zz, 


where YW (2/L") is another universal function which is related 


to B,(z/L") by 
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In a Similar fashion the temperature profile can be 


expressed as 
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Paulson (1970) has integrated (3-39) to get 
W, (2) | - wy) - 3 ie D, Pgs Cn [C1+6,)/2] 
(3-41) 
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Thus, having calculated 'z/L') from profile measurements, 


OD, (z/L"') can be determined by numerically solving (3-36) 


and YW, (2/L") can be obtained from (3-41). The velocity 
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44 
profile (3-38) is then specified for unstable conditions. 
The temperature profile can be obtained from (3-40) using: 
(3-20) and assuming that KL /Ky is constant. 

Another profile representation for unstable conditions 
was independently suggested by Businger (1966) and Dyer 
(unpublished). This scheme uses the interpolation formula 
(3-35) and the same hypothesis suggested by Pandolfo (1966) 
that (Z/L) = Ri. As was previously shown this leads to 
Deets @,(z/L). 


Now by definition of the terms involved 
‘ . 
Ales Syke (3-42) 


Substituting this into (3-36) yields 


) 
— 


Soyer’ byocesshe ape 
9) = (de viata 
Ginaset! dice yan! evel aed (3-43) 


From (3-20) and (3-32) 
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Businger-Dyer model and obtained 
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In free convection, as -z/L becomes large, the. 
Keyps representation predicts that both wind and temperature 
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vary asymptotically as z while the Businger-Dyer 
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representation implies V « Z and O0«¢z 


In the paper by Paulson (1970) it was found that 
the Businger-Dyer ae ntae ren gave the best fit to 
temperature profiles at Kerang,Australia and described the 
wind profiles for the same place equally as well as the 
Keyps model. There are more profile representations 


suggested in the literature but the two discussed here have 


received the most widespread acceptance. 
3.2.5 Actual Method of Profile Analysis 


It was originally intended that this thesis would 
compare the Businger-Dyer and Keyps profiles but a lack of 
data eliminated this possibility. The temperature profiles 


for August were completely lost because of equipment 
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malfunction. In October, the cold air over warm water 
created extremely unstable conditions, well into the 
free-convection regime. It was then possible to compare 
the observed profiles with the unstable limits of the two 
theoretical profiles but not with the interpolation models 
themselves. To this end a least squares fit was performed 
on temperature and wind versus the appropriate power of 
the height Z. 

In addition, flux estimates from the October 
free-convection profiles were desired for comparison with 
the flux estimates from the fast-response measurements. The 
heat flux was obtained using Priestly's free-convection 
temperature profile (3-27). A least squares fit was applied 
to temperature against z°-S which produced an equation of 


the form 
QO = Bet EO es (3-46) 


where both B and D are known. But comparing with (3-27) it 


can be seen that 
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(3-47) 


everything else in (3-47) is known and a profile estimate 


of the heat flux is available. 
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Knowing H permits Ky to be calculated using (3-28). 
Then, assuming that Ku = Kye the eddy diffusivity coefficient 


for water vapour, the flux of water vapour and hence of 


latent heat can be calculated from 


de (3-48) 


oz 


where E is the evaporative heat flux, e the absolute humidity 


= 


u 
|| as 
TL 

= 


and L the specific latent heat of evaporation. 

The shear stress € or friction velocity u* can be 
determined using either of the free-convection wind profiles 
Of Priestley ((oe30) or of Pandolfo*4{3-34). Considering (3-30), 
a least-squares fit of V and ares produces a constant F which 


must be defined by 
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u = i K M J H “4s (3-49) 
The term (elas y, 8,) can be taken from the temperature- 


profile equation which leaves K.,/K as the only unknown. This 


M 
ratio was already assumed constant in deriving (3-30). The 
value 1.3 from Lumley and Panofsky (1964) was chosen to 
complete the equation. 


A similar procedure fitting V to z - according to 


(3-34) yields another constant G such that 
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where C, = 3k C and C is the same constant as the C in 
(3-47) for which there are values already calculated. This 
method does not require the dubious assumption of a constant 
Ky/Ky- 

All derivatives were approximated by 
26) ae F ur ) tn (2) 
oz (2,2, )% (2, 2) * a, 
following a suggestion by Panofsky (1965), where F is a profile 
variable. The approximation applies at the geometric mean 
height of the observation levels and is rigorous for 
logatathmic! profiles’ (i.e. neutraljconditions)/: For Fez? 
the error was estimated by Paulson (1970) to be ~6% for 


observation heightsof 0.5 and 16 m.and should be less for 


the heights used in the present study. 
3.3 “Vertical Fluxes by Eddy Correlation 
Shesd alntroduction 


It is a standard procedure when examining the 
characteristics of turbu-ence and turbulent transport to 
assume that all quantities subject to variation because of 


turbulence can be written as the sum of a time mean anda 
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fluctuating part (see e.g. Lumley —_— Panofsky 1964), 
where the mean of the fluctuating part is zero. In the 
present study an overbar will indicate mean or average 
values and the prime will indicate the perturbation from 
the mean which has zero average. Thus the instantaneous 
value of some physical parameter X ina turbulent flow can 
be expressed as 
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SOT om or ne OX, i350) 


2,a6e ‘Turbulent. Fluxes 


The equations governing flow in the surface 


boundary layer are given by (see e.g. Lumley and Panofsky 


1964): 
a 
du Uy PL Pee a de yy ou: 
dt d x; Po Ok ats he 
ee ae 
dx; Vero) 


2. 
PR 5 PK a ) ane 
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dt aX; oy BMetoC? 
These are the equations of motion for a viscous, Newtonian 
medium in a uniform gravitational field, neglecting rotation. 
Subscripts imply summation over the three coordinate directions 
and x, is the vertical axis. The fluid speed is denoted by 
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viscosity and is the thermal diffusivity. The temperature 
T is a deviation ee the name rarned adiabatic value denoted 
by T,. In the same way P refers to a small deviation from 
the undisturbed atmospheric pressure PR, and lai refers toa 
deviation from the undisturbed atmospheric density fo ° 


The Kronecker delta, iG 


ee is unity if the indices are equal 
but otherwise is zero. Motions at very low Mach numbers in 

an atmosphere consisting of a perfect gas of constant 
composition have been assumed, and variations in the viscosity 
and thermal diffusivity were neglected. 


Writing the velocity, temperature and pressure in the 


torm of (S-oL) 


T = T 2 a (F353) 
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and substituting into (3-52) yields: 
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ag" JE dx; Ox; dx; dx, 
jad, re. OP gy dU; + py? u; Pa i Kat &, 
ee ps I; dx aba ap : 
ou, 4 du =, 0 
dj dx; (3-54) 
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Taking averages produces the equations for the mean motion 


of the turbulent fluid: 


De om tl ome dP: oo, 13 T Ss 
dt I; dx; See pacer ey: 
OU: Nhe 2, 
IX 
(3-55) 
JT 1 UzdT 4 U; > Doh = Jd Pumie 
t dx; o xe P) x; 
Subtracting the middle members of (3-54) and (3-55) 


gives du, Hon, 70 and the term u; P) u; VARS 


to Wx; ( ae u;” ) 


can be transformed 
Multiplying this by the density fe 
produces the divergence of a momentum flux which represents 


a force due to the turbulent transport. of momentum... This 


term is frequently written on the right-hand side of the first 
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of (3-55) and added to the stress giving 


Ch UV dU; ee. gree (3-56) 
dt dN; Ving Ox, iy 


Te ovat tofu? pale +25 | yarittey 4 (3-57) 
Ix, 

The analogous term in the temperature equation 
u: dT, On; when multiplied by Cp Po can be transformed in 
the same way to give the divergence of a heat flux, 
representing the heat transferred by turbulent temperature 
fluctuations. Similarly this term can be taken to the 
right-hand side of the third member of (3-55) and a heat 


flux vector defined as 


~ —_ / 
2 ais a C a oa 
Hy =~ ¢p }{ 97 + opps (3-58) 
dX: 
p) 
Neglecting molecular viscosity and diffusivity and 


assuming no external pressure gradient (3-57) and (3-58) 
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Hj = Cp /o u; 


Finally the assumption of horizontal homogeneity 


reduces the expressions for the shear stress and heat 
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Pix to 


sy fas Pe Ler (3-60) 


Hity sam Beoget wie Eb (361) 


Here u' is the perturbation of the horizontal wind and w' 
is the perturbation of the vertical velocity. 

The vertical turbulent flux of a quantity then is 
determined by the covariance of the perturbations in the 
quantity and perturbations in the vertical wind. If, for 
example, upward fluctuations in the vertical wind correspond 
to an increase in X or a positive X' and downward fluctuations 
occur with a decrease in X then a net upward flux of X results. 
The easiest way to visualize the significance of the 
covariance is to consider that upward vertical motion is 
bringing air with a certain X value up from somewhere below 
and downward vertical motion brings air with another value 
of X down from above. If there is a gradient of X with 
height then higher values of X will be carried to levels 
having a lower value of X and a flux of X "down" the gradient 
or xX results. Inen, corresponding to (3-60) and’ (3-61), if 
the vertical wind is denoted by w= w+ w', the turbulent 


flux of X is given by 


egies x" (3-62) 
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The perturbations (e.g. w' and X') were not measured 
directly for the evaluation of Aes eat of the form of 
(3-62). However, the covariance of the perturbations, and 
hence the flux, can be related to the quantities that were 
measured, namely the total quantity (w and X) and the mean 
quantity (w and X). Using (3-51) the covariance of the 


perturbations can be related to the measured quantities by 


eX Be ete ar eae ae 
Seine un) Aueesaneity ain nd 
See te aN, on ee wT 
= veRuGe wand wie ks (3-63) 


In stationary,homogeneous flow w = 0 and the second 
term on the right is zero. The turbulent fluxes of heat 
(H = eee water vapour (E = w'e'), and momentum 
(-ur? = u'w') can be calculated if values of vertical wind 
w, temperature T, water vapour content e, and horizontal wind 
velocity u are measured with sufficiently fast-response 


instrumentation to make the covariances with w accurate. 


3.3.3 Fluxes Associated with a Mean Vertical Wind 


— eS Eee 


If the mean vertical velocity is not zero then the 
conditions of stationarity and horizontal homogeneity no 


longer hold and the turbulent flux at a certain height cannot 
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be considered representative of the flux from the surface. 
The fluxes due to the mean vertical wind can be estimated, 
however, and a better estimate of the surface flux can be 
obtained by adding the turbulent and mean wind fluxes. 

The contribution to the vertical fluxes by the mean 
vertical wind arises from the second terms of the first and 
third equations of (3-55). The term U; Bey, Joe can be 
taken to the right hand side of (3-56) and incorporated into 
the shear stress YT using the equation of continuity in 
(3-55). Assuming horizontal homogeneity and neglecting 
molecular viscosity and horizontal pressure gradients, the 
expression for the shear stress, including the contribution 


from the mean vertical wind, becomes 
i fie: LEPAP Ae: 34 fo ww (3-64) 


Similarly, the expression for the heat flux, 
including the contribution of a mean vertical wind, is given 


by 


=) 


H = Cp Pe Wo lee ea Cp Po W (3-65) 


where again T is the deviation from adiabatic conditions 


caused by the heat added at the surface. 
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CHAPTER TV 
RESULTS 
Au. LNeroduction 


Before any measurements could be properly interpreted 
the reliability of each had to be estimated. This was done 
using available literature on the instruments and calibrations 
carried out at the University of Alberta. Tables 4 and 5 
give the estimated errors of measurement of the slow-response 
and fast-response instruments, respectively. Errors caused 
by terrain slope do not apply here of course because the 
lake surface is always nearly horizontal but errors in 
aligning the fast-response instruments with the vertical 
must be considered to be of possible significance. Such 
errors were estimated to be less than 15 minutes of arc for 
the transit procedure used at Lake Wabamun. 

If @ is the deviation from the vertical then the 
error in vertical velocity is approximately u sin9 or 1=2 
cm/sec for the horizontal wind speeds of interest to this 
project. This does not appreciably affect the turbulent 
fluxes because the mean vertical wind is subtracted from the 
flux estimate. However, estimates of the mean vertical 


velocity and associated fluxes are definitely influenced. 
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Sensor 


Cup Anemometer 


Thermistors 


Hygrometers 


Thermistor 


TABLE 4 


Slow-Response Instrument Error 


Element 


Wind Speed 


Air Temperature 


Dew-Point 


Water Temperature 
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TABLE 5 


Fast-Response Instrument Error 


Sensor Element Eevor 
Sonic Anemometer Vertical Wind Speed agt 
Propeller Anemometer Vertical Wind Speed eo 
Propeller Anemometer Horizontal Wind Speed 10% 
V component as much as 
; 3037 
Sonic Thermometer Air Temperature 53+ 
Resistance Thermometer Air Temperature 53 
Lyman - « Humidiometer Vapour Pressure 53 
Sonic Anemometer 
Thermometer Heat Flux 103+ 
Fluxatron Heat Flux 1 5 ig 253° 
Sonic-Lyman - « Vapour Flux g 
Shear-Stress Meter Momentum Flux 3° 
L Kaimal and Businger (1963). 
2 B. B. Hicks (1972): For -wind-speeds of interest’ here. 


i McBean (1972): 
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Largest error for near-neutral conditions. 
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29 
The shear-stress meter consisting of three orthogonal 

propeller anemometers was built so that the axes of the 
mount were within a few minutes of arc of being mutually 
orthogonal. Any error from non-orthogonal measurements was 
considered insignificant compared with the error introduced 
by the poor response of the propellers at large angles of 
attack to the wind. This latter error was included in the 


values given in Table 5. 
4.2 August Tests 
4.2.1 Local Weather Conditions 


As was previously noted the weather conditions for ' 
the two days of data collection in August, 1972 were quite 
Similar except for wind direction. The Edmonton International 
Airport, fOr ky miles east of the lake, recorded maximum 
temperatures of 80 F and 79 F on August 24 and 25, respectively. 
Wind speeds for the two days averaged 3 m/sec and 2 m/sec 
with the greater speed occuring on August 24, the day of 
southeast winds and measurement of fluxes over the heated- 
water plume (see Figures 4 Ana 5). Scattered cumulus clouds 


were observed on both days. 
a. 2.2 stability 


Meaningful Richardson numbers could not be calculated 
for the test days in August because of the unavailability of 


temperature profiles. However Z/L was calculated for 
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FIGURE 6, Surface temperature upwind from tower site; Lake 
Wabamun, October 21, 
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the two days using equation 3-2. These data, presented in 
Table 6, indicated near-neutral eendi ions in the flow 
across the unheated lake and instability on the day with 
flow across the heated plume. The Z/L values were based 
upon thirty~minute means of u* from the shear-stress meter 
and heat fluxes for the same intervals from the sonic 
anomometer-thermometer, hence the valid height is 3 nm. 

The values of u* are suspect for August 24 (see 
Table 7) because the vertical propeller of the shear~stress 
meter was not recording. The only record of vertical wind 
was from the sonic anemometer which was separated in space 
from the shear-stress meter and has a different response 
time. An average value of u* was used to calculate all of 
the L values in order to minimize the effect of these errors. 

Typical wind profiles, averaged over 30 minutes, are 


presented in Figure 7. 
4.223¢9F Luxes 


The fluxes for the two test days in August, calculated 
from the fast-response instrument data, appear in Table 7. 
Over the plume (August 24),the heat fluxes were on the average 
considerably larger than those over the rest of the lake 
(August 25), while evaporation from the plume, as measured 
by the Lyman-alpha humidiometer and sonic anemometer, was 


3 to 5 times greater. The inconsistent momentum flux on 
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TABLE 6 


August Z/L Values 


Date/Time Wind Speed Sonic 2/L Flax’Z/L 
(cm/sec) 

2471515 452 -0.0068 
1530 386 ~0.0101 
1545 397 -0.0096 
1600 378 -0.0170 
16235 353 -0.0892 
1630 341 -0.120 

256 £330 £O3 =0.0020 -0.0011 
1345 149 -0.0080 -0.0010 
1400 “4 136 -0.0072 -0.0013 
1415 144 -0.0022 -0.0026 
1430 160 -0.0013 -0.0015 


1445 £99 -0.0034 -0.0006 
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67 
August 24 was explained in 4.2.2. 

The difference in fluxes on ne two days is probably 
not entirely accounted for by the heated-water plume. Part 
of the decrease in magnitude of the fluxes on the second day 
can be explained by the lower wind speeds observed on that 
day. The lower wind speed would indicate less mechanical 
mixing and smaller fluxes. The difference in surface 
temperature was undoubtedly the dominant factor, however. 

For August 25, with heat flux estimates from both 
the sonic equipment and the fluxatron available, a marked 
discrepancy between the two heat eta measurements existed 
as can be seen in Table 7. In almost every case the sonic 
estimates were higher than those of the fluxatron. One 
possible explanation is linked with the fact that the response 
time of the sonic anemometer is considerably less than that 
of the propeller of the fluxatron. This would result in an 
underestimate of the heat flux by the fluxatron. However, 
this should account for only 25-30% (See McBean, 1972) of 
the ob ears difference and a satisfactory explanation of 


the total discrepancy is unavailable. 
4.3 October Tests 


Joel 1Stabdhity 


On October 21, 1972. there was a flow of cold air 


(+1°C) which crossed the heated plume before arriving at 
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the instrumentation site. Unstable conditions were 
expected in such circumstances and ae Richardson number 

and Z/L' as calculated from the profile measurements 
indicated extreme instability. Indeed, as Table 8 shows, 
[Ril as calculated for 1.3 m was greater than 1.0 throughout 
the data collection period and the surface boundary layer 
was therefore well into the free-convection regime as 


aefined an) 3.223. 
A.3a2 October Profiles 


Because of the definite free-convection situation, 
an attempt was made to fit the theoretical free-convection 
profiles (see 3.2.3) of wind and temperature to the data. 
The fit was made for half-hourly mean values and the results 
are shown fin Tables 9 and 10 and in Figures 8 and 9. Figure 
8 shows observed temperature (crosses), the 1/3 power law 
(solid line), and the 1/2 power law (dashed line). The 1/3 
law, of course, was predicted by the Priestley free-convection 
theory and is the value to which the Keyps formula is 
asymptotic. The Businger-Dyer theory is asymptotic to the 
1/2 power curve. Both curves fit the data well but as is 
shown by the sums of squared differences in Table 9 the 1/3 
power, law was slightly better in every case. 

The fit of the theoretical wind profiles to the 


measured was not as good because of a larger uncertainty 
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FIGURE 8, October temperature profiles, Observations are 
indicated by triangles,The solid curves represent the 173 
power law aud the dashed curves the 1/2 power law. 
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(see Table 4) in the wind measurements (see Figure 9). 

Once again, as is shown in Table ihtetna departures of 
theoretical values from actual were all within the error 

of measurement although thé Pandolfo 1/6 power law seemed 

to fit the observed values to a slightly better degree. The 
1/3 power law corresponds to the Keyps free-convection 
prediction again while the Businger-Dyer formula is asymptotic 
to the 1/4 power law. 

Taylor (1970) described a numerical model of a flow 
over a step change in surface temperature, a situation 
analogous to the actual conditions neeiene Wabamun on 
October 21. Comparison of the observed profiles with Taylor's 
results was difficult because Taylor presented all predictions 
in terms of deviations from an upstream neutral profile. 

At Lake Wabamun, no information on the profiles upstream 

of the heated water was collected. The numerical results 
indicated an acceleration of the flow in the lower levels to 
values above those at the same height upstream, with the 

fluid peehae wae ae propagating upward with increasing 
distance downstream from the temperature change. The observed 
winds cat 0.5 m. did indicate a trend in this direction as 
shown in Figure 9 but the deviations from the predicted 
profile were within the expected error of measurement of 

the anemometers. In addition, for the actual distance 


downstream of the increase in surface water temperature 
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@ 
(300 m), Taylor's model predicted profiles adjusted to the 
new surface to a height above the chose level on the tower. 
In short, the whole profile should have adjusted to the new 
surface conditions. It is possible that, if the low-level 
acceleration was real, it was the result of the air flow 
meeting higher temperatures as it approached the tower site. 
The warmest water was along the shore as shown in Figure 6. 
Further comparison with Taylor's numerical results would 
xeally only be speculation, 

Calculations of the fluxes of momentum, heat and 
water vapour from the profiles were carried out as described 
in 3.2.5. The derived values of u*, and of sensible- and 
latent-heat fluxes appear in Table 1l. The dew-point sensor 
results were unrealistic because of the dense fog conditions. 
The latent ‘heat fluxes were obtained by assuming saturated 
air up to 7.5 m, equating air temperature to dew~-point 
temperature and substituting into equation 3-48. This 
assumption seemed justified by the presence of fog extending 
above the eae level of measurement and by the observation 
of condensation on instruments and other equipment at all 


levels. 
Ay 333 Directly-Measured Fluxes 


The turbulent heat flux as measured by the sonic 


anemometer-thermometer and by the thermistor and propeller 
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Momentum 
Flux (mw/om~) 


—— 


Pandolfo Profile 


October 21 

30 Mins Sensible Heat Latent Heat 
aerate Flux (m¥/om) Flux (mw /em*) 
2009 8.40 8.88 

2015 S%a2 8.64 

2030 3.05 BZ 

2045 1426 G37 

2100 RAS 6.87 

ZES 7.47 7.30 

2130 7? ou y a ae 
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anemometer were in good agreement with each other and with 
the profile estimates for the Bi Perens running means 
shown in Figure 10. There was a gap in the fast-response 
instrumentation data and the discrepancies in the last 
hour's results between fast-response and profile estimates 
were based on a rather large heat flux in the final fifteen 
minutes as calculated from the fast response data and not 
duplicated by the profile estimate. Generally the fast- 
response heat fluxes were within 10% of each other and 
within 15% of the profile estimates of heat fluxes except 
for the final 30-minute period. 

As can be seen in Figure 1l the latent-heat fluxes 
as measured by the Lyman-alpha humidiometer and sonic 
anemometer were generally 30% lower than the values obtained 
from the profiles. There are at least two possible reasons 
for the difference. The humidity profile was based upon 
the temperature profile assuming saturation, but despite 
the fog this may have been an incorrect assumption. Also, 
the profile water-vapour flux was calculated assuming that 
the eddy coefficients for heat and water vapour transport 
were equivalent which may not have been the case. The vapour 
content as measured by the Lyman-alpha would produce a 


relative humidity of onl 60%/which is ajseriousbno as yet 


unresolved inconsistency. 


need VER, Fedtd odd oak! syPS oad opted atoeiyaix” e. 


xVoqRY ante | ae ae new ove ‘ton ye bait ur tneisviepe all a he 


ey | ue we 


sis Ww iris TON fone “aoe 


arkoyage _ er nd ‘wag a ae aeot | OL 4 | 
aul ae ae 

ee ae Y eoionnqnoelt it bie aieb nolan 

ao sariaas ot iteng bas abtssianiaasiil nvemgod 


Say: Siege Este DeneqRer tae2) ants a, osmlamne ea 
~seu> off Yi exeaen” oremiree olziiva adda yal bs 

hai <onise ions 26 #02 niki we wae sec Ta isi es ms 
fqame eoxut® peat to antamiaeg saad ons Ro wet 


eee Bets snaand. ort a oni nk aoge ad mad eh wh 
elt: bite’ 1 omg. tiniot grea nite ws 
banttintde sallow faa tists ‘mewor BOE <a " | 
siieidi oiiteeem, oa aesiot, we axe ecb | 
a. boned, waste Ecsicud eibeat — 
Naknie isso eit ae rived eared’ ames por may, * 
40h) Bi TRE hyo kacitee Baw xi, iuoqev~ xed 6M etttoxg ofa 
saoqaner? parca “fodew be sort 4@2 ajnoisi*ises: YbbS. ota 

~g toubose dies tg tng edt ya bosueeisay ae Saedgd 

Joy8b Bi amp tiene! wt! Abi MOD “tno Yo yrthinad owkgeter 


4a 
. ; e 
aphodatanten Bovicneznd 
bs Pea ann) crs ae if 
oF iad « ie Re 5 ; ne 7 + 5 
- i : + . f 
Fe i pat \ s “s iY 
a Wotan 
. fi oaths eg 
2 : : | sa PANT. ha 
‘ 421 Pony ran i; 
; 1 ay hia sah 
| 5 ED ee ee 
“ae i mye Bt cr MN Pas any, 


80 


aie squtod peqjotd [ly °o Aq elep ITUOS 
peqyeoTput o1e ejep oeTtTjoad worzz poeyeTNoTeo 


*sSUuBPOW DINUTU-(E UO posed 
woij osoy} pue x Aq e2ep VoAeXNT} Worzzy Ysoy ‘oe Ag 
SOXNT A “19q0700 ‘suostieduod xnTy eoYH ‘OL HeNDIa 


QWLL LHOLLAVG NIVINNOW 
oOl2 O£O03S C002 O€6l — O06 
= ik a 71S 
=) 
° 2 


\ 
\ 
\ 
\ 
\ 


N 


we) 


= 


NHS 


is 


XO1TI LVGH ATL 


(mo (a 


06 


ry Le wali a i iid hal 


gq depres 2c5. 


a2 ph gai byosts 


* 
ae 


Rott. q¥cs gh s oy ppaN4 {10m sauge @ 


f COMES LTP. 2 


' 
ica 
ae 
4 
= 


*. GSEpveh” Eyrese asyperrered grow Sperpye qyes exe. poqrescéq 


hau te lial. dy aut ibaa 


ES 


~ 
~ 


be 


\ 


<6 weieag® 


Ei i 


is a # . wire A Vr : ee 7 } 
i - he > : fi t & &. A . ee ae nie 4 
ais : ' J ¥ ; , i 7 i ; , ‘i 
: oe ue ee 1 A ne IRR 
Reh a. | aie ste Oy Y als r i : } ; 


a ec, Oe nn eens Hc 


* 


FG 


pees: 
a 
bEGISE 


= 
4 
La 


i if ce 
a 


: : > ; AFD ¢ 7 a j PN: i Pin 

7 . 3 a ee i a iy ier APA ai). ; a 
ih F it wi Sie r i : P = 4 a oe at at i] 
iw § i ~ ps | ‘ if fi 7 

= : : ee tings | i ri ys ee Pe a | - - 


81 


aie squtod poeajotd 
aie eqep eTtjoad wory 


*supow oy NnuTw-OF 
tty °o Sq eqep Jueunazjsut osuodsei-jseyJ wory Eso pur x Aq 
poqyetnoyeo sexnzq “leqoj909 ‘suostiedumoo xnTZ Ie9Y JUuSIeT 


GALL LHOLIAVA NIVINNOW 


uo poeseq 
pezeoTput 


“TT geno 


29) 


~ 


© 


co) 


(_wo/ma) XOId LVEH LNALVt 


c 


i 


Sehorne * 


+= 


~~ 


y 


Reb x82 


; ' 
' 
3 
4 
he, 
* 
i i 
) ee 
we 
Ce, 
ve 
Penny 
4 f 
fi 
Ps | a : 
eT | . Fh sbhal oy 
; ia 
* ’ ey 7 
i A ) = 
| U) : $ ‘ - g = 
‘ 14 
' J rn 
me ) | 4 a * os ea! vt 
: 7 : 5 r ae Pe Gs) 
ae i 
’ + a | ' Lae | ig 
* “ > G ix is ae 
ideas ho ne H > al gee u's ae 


— 
*~ 
6 

72 


Cle 
aeeiies et 
We 


DeAL 
Cy 
. 
‘a a 
: 
: 
-_— 
—s = 
i: 
. 
ost 
os 
= 


i Vanes * i 


eee i. ee | 
ee in uy \. 


ie 
VE 
os 
es . 


74 ‘node ra, 


eo 4 ; ai 
7 ny 4 


fi 


82 
The momentum fluxes as calculated by the shear-stress 
meter were even more surprising. Not only did the shear 
stress u* differ markedly from the profile values of u* but, 
as can be seen in Figure 13 the momentum flux, or -u*4, was 
frequently upward against the velocity gradient. In fact 
the average momentum flux over the entire data-collection 


period was upward. 
4.3.4 The Anomalous Momentum Fluxes 


Figure 13 shows five-minute averages of ~ur? Or Ot, 
the momentum flux u'w' as measured by sthe three-comoonent 
shear-stress meter. Momentum fluxes calculated from averages 
over periods as long as forty-five minutes did not show a 
difference in sign from that shown in Figure 13. The period 
from 115 to:-150 minutes (2025-2100 MDT) was plotted as a 
horizontal line because no data were available for that 
period. ~ Mniorder jto further investigate the’ cause of: the 
upward momentum fluxes a few five-minute periods were chosen 
forma closer examination. The two largest upward spikes 
corresponding to 10-15 minutes (1840-1845 MDT) and 50-55 
minutes (1920-1925 MDT) were considered to be the obvious 
Eaten suitable for more intensive study, while the 
downward fluxes in the five-minute periods preceding each 
positive spike were used for comparison purposes. In addition 


several other intervals throughout the data collection period 
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FIGURE 12, Five-minute heat fluxes (sonic); 1830-2130 MDT, 
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84 
were chosen randomly to be examined more closely. 

In each chosen interval, means of total horizontal 
wind u and of vertical wind w from the shear-stress meter 
were examined as well as the fluxatron temperature T. The 
original intention was to find a time period over which to 
average, large enough to eliminate the fine-scale fluctuations 
but short enough to isolate the correlation of u and w 
necessary for the upward momentum flux. Normally of course, 
the horizontal wind increases with downward vertical velocity 
as faster moving aix from above is brought down and vice 
versa for upward vertical velocity. For an upward momentum 
flux the opposite must occur at some scale and upward vertical 
velocity must correlate with an increase in horizontal wind. 
After some trial and error, thirty-second means were found 
to show best the anomalous correlation between vertical and 
horizontal velocities. 

Figures 14 £0°20 show the plotted values! of the 
thirty-second means. In each case u is the top curve 
(stars), w is the bottom curve (squares) and temperature T 
(crosses) is plotted in the middle. The vertical scale 
represents cm/sec for u, (cm/sec)/10 for w and deviation of 
temperature in °K/100 for T. Every test interval showed both 
the normal situation with u increasing with downward w and 
the anomalous situation of u increasing with upward w. Some 


of the more accentuated examples of the latter are indicated 
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FIGURE 14, October 30-second means; 1833-1845 MDT, 


| S| 
(a9) 
=) pp pt pp tp 
, | 
: | 
rH 
S = x ' 
5 \ + 
: ye. | 
. t 
2. 
cit ue 
cS 
a 
=D) 
9} — 9 — 
oO : 
=| | 
U.00 3500 6.00 9.00 12.00 15.00 


TIME IN MINUTES 


FICURE 15, October 30-second means; 1900-1913 MDT, 
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FIGURE 16, October 30-second means; 1916-1924 MDT, 
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FIGURE 18, October 30-second means; 1938-1950 MDT, 
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FIGURE 20, October 30-second means; 2118-2130 MDT, 
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89 
by vertical lines on the figures. All intervals show the 
expected correlation between w and T. Using these thirty- 
second means, suitable one-minute intervals were selected 
for which the fine~scale fluctuations were examined. 

The one-minute intervals thus chosen were parts of 
the five-minute periods preceding and making up the two 
largest upward momentum flux spikes in Figure 13. For each 
minute chosen, half-second means of all the instrument 
output were computed and these values are plotted in Figures 
2. through 49. The total Uis.,the, total horizontal wind 
computed from the horizontal Lae ears of the shear-stress 
meter. This in turn was broken down into a one-minute mean 
direction component and a component perpendicular to the 
mean direction, the latter being denoted by V. It should be 
noted that anergy Gave and sonic vertical velocity values 
are given relative to some arbitrary zero. 

The two minutes 1835-1836 and 1919-1920 (Figures 21 
to 25 and 38 to 42) are parts of five-minute periods which 
yielded the expected downward momentum flux. In both minutes 
there was a strong Legere pe between increasing upward 
velocity and decreasing horizontal wind. As well, both show 
a strong positive correlation between the sonic temperature 
and vertical velocity indicating a large positive heat flux. 

The period 1842-1843 (Figures 26 to 31) falling 


within the first large upward momentum spike, was probably 
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FIGURE 22, October half-second values; 1835-1836 MDT: sonic 
temperature, 


oO 
r=) 
oO 
a a a 
| | 
= | 
ca i 
+. 
2 | 
Li 
W9 
eo | 
LIS ‘A | 
ag ee _ iS ‘ia 
2 Svan | rw \ ee 
r——4 c¢ aw 
cis al as 
oO 5 
Mo 
= wi 
0 fp pt tt tt tt 
0,00 12.00 24.00 35.00 48,00 60,00 


TIME IN SECONDS 


FIGURE 23, October half-second values; 1835-1836 MDT:sonic 
vertical wind, 
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’ FIGURE 24, October half-second values; 1835-1856 MDT:total 
horizontal wind, 
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FIGURE 25, October half-second values; 1835-1836 MDT: shear- 

stress meter vertical wind, 
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FIGURE 29, October half-second values; 1842-1843 MDT: shear- 
stress vertical wiad, 
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FIGURE 31. October half-second values; 1842-1843 MDT: horizontal 
wind, mean direction component, 
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FLGURE 32, October half~second values; 1843-1844 MDT: sonic 
temperature, 
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FIGURE 33, October half-second values; 1843-1844 MDT: sonic 
vertical wind, 
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FIGURE 34, October half-second values; 1843-1844 MDT: total 
horizontal wind, 
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FIGURE 35. October half-second values; 1843-1844 MDT: shear- 

stress vertical wind, 
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FIGURE 30, October half-second values; 1843-1844 MDT: horizontal 
wind, V component, 
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FIGURE 37. October half-second values; 1843-1844 MDT: horizontal 
wind, mean direction component, 
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FIGURE 38, October half-second values; 1919-1920 MDT: total 
horizontal wind, 
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FIGURE 39, October half-second values; 1919-1920 MDT: shear- 
stress vertical wind, 
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FIGURE 40, Octeber half-second values; 1919-1920 MDT: horizontal 
wind, V component, 
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FICURE 41, October half-second values; 1919-1920 MDT: horizontal 
wind, mean direction component, 
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temperature, 
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FIGURE 44, October half-second values; ;920-1921 MDT: sonic 
temperature, 
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FIGURE 45, October half~second values; 1920-1921 MDT: sonic 
vertical wind, 
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FIGURE 46, October half-second values; 1920-1921 MDT: total 
horizontal wind, 


100.06 


SHEAR W CM/SEC 
0.0c 
? 
4 


un oer ee 
'o.00 12.00 ou, 00 36,00 U8, 00 60.00 


TIME IN SECONDS 


FIGURE 47. October half-second values; 1920-1921 MDT: shear- 
stress vertical wind, 
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FIGURE 48, October half~second values; 1920-1921 MDT: horizontal 
wind, V component, 
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FIGURE 49, October half-second values; 1920-1921 MDT: horizontal 
wind, mean direction component, 
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the most interesting of all. The curves of shear stress w 
and horizontal wind u were strikingly similar and in phase, 
each having two fairly broad peaks within the 60- second 
micerval “(See Figures 28 ahd 29). Thus-on this large “scale 
the horizontal velocity increased as upward velocity increased. 
A close examination showed evidence of a finer scale 
fluctuation, superimposed upon the 30- second period 
oscillation described above, which had the expected increase 
of horizontal wind with increasing downward vertical velocity. 
Considering the larger scale again, the V component in 
Figure 30 showed a definite change in flow direction at the 
same time as maximum upward velocity was achieved. The sonic 
temperature and vertical wind were in phase for the first 
30 seconds but were almost perfectly out of phase for much 
of the second 30 seconds. 

The next minute, 1843-1844 (Figures 32 to 37), showed 
virtually none of the large-scale oscillation. Only the 
fine-scale fluctuations corresponding to the usual downward 
momentum flux were apparent. The sonic temperature and 
vertical wind traces showed no strong positive or negative 
correlation. 

The final minute, 1920-1921 (Figures 43 to 49), was 
also taken from a period of large upward momentum fluxes and 
the corresponding data were similar in appearance to those 


for the interval 1842 -to 1843, although the large-scale 
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fluctuations were not as easily separated form the small 
scale. Reduced mean wind speed may have been an important 
factor in this regard. The V component indicated no truly 
Significant direction change (Figure 48) and the sonic 
temperature and vertical wind were in phase (Figures 44 and 
45). 

In brief then there appeared to be a large-scale 
turbulence somewhat random in occurrence at the measuring 
Site that accounted for a large upward transport of momentum. 
Superimposed on this was a small-scale turbulence continuously 
operating and transporting momentum in the conventional sense, 
When only the fine~scale turbulence occurred there was a 
small negative momentum flux but when the large-scale and 
small-scale turbulence occurred simultaneously there was a 
net upward transport of momentum. One can speculate that 
the large-scale phenomenon was some sort of organized 
convective eddy tilted in the flow direction to enable it to 
transport momentum against the mean velocity gradient whereas 
the fine-scale fluctuations represented mechanical turbulence. 
There was no evidence of a change in the temperature during 
the passage of the convective eddy and only a suggestion of 
Some sort of turning of the horizontal wind, a rotation, 
during its.occurrence. 

It should be noted that Figure 12 shows strong 


downward momentum fluxes in the conventicnal sense for the 
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period after 150 minutes (2100 MDT). Higher mean wind 
speeds in that period may have been responsible. It is 
conceivable that the stronger winds destroyed the organized 
convective cells, leaving a small-scale, more nearly 
isotropic turbulence accompanied by downward momentum 
transport. 

Others have found similar counter-gradient momentum 
fluxes during highly unstable conditions. In particular 
Kaimal and Businger (1970) described in considerable detail 
such an occurrence in Kansas. They examined an irrotational 
convective plume and a rotational dust devil, both of which 
accounted for an upward momentum flux. The analogy between 
the Kansas and Wabamun data seems to end here, however, 
because the upward momentum flux in Kansas was accounted for 
by a burst of upward flux in the short-period or fine~scale 
turbulence. The Wabamun data on the other hand showed 
consistent downward fluxes by fine-scale turbulence. An 
upward flux was evident only when periods longer than 30 
seconds were considered. Half-second and one-second average 
fluxes for the Kansas occurrences showed considerable upward 
momentum flux but one-second average fluxes for the Wabamun 
data in the time period 1842-1843 showed all downward fluxes. 

The Kansas data were taken at higher levels (5666 om 
and 22.6 m) but this should only have brought them into a 


region of Richardson numbers similar to those experienced at 
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lower levels in the more unstable conditions existing at 
Lake Wabamun in October. The other difference of course was 
that the convective systems moved past the observing site 
in i aepeicteiit at Lake Wabamun the site was fixed to the 
source of the convection, the hot water plume. An extensive 
observing set up enabled Kaimal and Businger to deduce a 
description of the Kansas convective systems but the limited 
data at Lake Wabamun prevented similar deductions about the 
structure of the Wabamun convection. The only possible 
conclusion is that the Kansas and Wabamun convective systems 


were dissimilar. 
“23.5 Turbulent Flux: of- Liquid Water 


Because of the high concentration of liquid water 
(fog) in the air during the October experimentation an effort 
was made to estimate the vertical flux of this quantity. 
Because the liquid water had evaporated from the plume 
surface and then recondensed it represented a legitimate 
Ebux tof Te wank heat from the surface. 

No direct measurements of this quantity were made 
and several assumptions were necessary. The only indication 
of the quantity of liquid water present was based upon the 
visibility which was approximately 300 £t. Wiener et alk, 
(1961) quoted a liquid water content for this visibility in 


ocean fog of 0.3 g/m? and this figure was accepted for the 
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Wabamun situation. This was an order of magnitude smaller 
than the water-vapour content ree ee at approximately 
“3 g/m? by the Lyman-alpha humidiometer. Then, assuming that 
the gradient of liquid water was the same as that for water 
vapour and the transporting mechanisms were similar, the 
flux of liquid water was also an order of magnitude less 
than that of water vapour. This would account for a 10% 
increase in the total flux of latent heat but because of the 


uncertainties involved in calculating the liquid water flux 


no further attempt was made to include it in the results. 
4.3.6 Fluxes Associated With the Mean Vertical Wind 


The mean vertical velocity at Lake Wabamun in 
October was not zero and in fact an upward velocity 
approximately six times larger than the possible instrument 
error was measured thorughout the data collection period. 
As previously discussed, the existence of a mean vertical 
motion indicated that the turbulent fluxes no longer represent 
the total flux from the surface since the mean vertical wind 
carries a significant part of this total flux. 

The vertical transport of quantities resulting from 
this mean upward motion were calculated as outlinedin 3.3.3. 
The temperature of the undisturbed atmosphere was chosen as 
the temperature at 7.5 m for the purpose of calculating the 


temperature deviation from adiabatic because the temperature 
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profiles seemed to be approaching the adiabatic lapse rate 
at this level. If this temperature deviation is denoted 


by T then the flux of sensible heat from the surface due to 


the mean W is given by (See 3.3.3) 


Ww 


H = ¢ p Tw 


The vertical flux of water vapour associated with 
an average upward wind is the product of water vapour amount 
and vertical wind at that point. The situation is analogous 
to the heat case, however, and not all of this is a flux 
from the surface. Most of the moisture is brought into the 
area of upward velocity from elsewhere by convergent winds. 
The water vapour at a level above the heated water is therefore 
composed of what is advected in horizontally and what is 
added from'below. The flux from the surface associated 
with the mean vertical motion is the product of the vapour 
added to the air by the surface and the mean w. The amount 
of vapour added by the surface was handied in a manner 
similar ene way the temperature deviation was calculated. 
As previously mentioned, the temperature of the undisturbed 
air away from the plume was assumed to be that of the 
7.5 m level. For lack of information and to provide a 
minimum flux estimate, this air was assumed to be saturated. 
The estimated amount of water vapour in the undisturbed air 


was subtracted from the amount of water vapour in the air 
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at the tower site to produce an estimate of the water vapour 
added by evaporation. The vapour Fax from the surface 
associated with the mean vertical wind was then calculated 
by taking the product of this water vapour difference and 
the mean vertical wind as measured by the shear-stress meter. 

The fluxes claculated in this way, along with the 
turbulent and total fluxes of sensible and latent heat are 
presented in Tables 12 and 13. In each period the sensible 
heat amounts transported by the two methods were of similar 
magnitude while the iatent heat flux because of the mean 
vertical velocity was considerably less than its turbulent 
counterpart in every case. The rather large uncertainty in 
the vertical velocity (~2 cm/sec, see 4.1) produces a 
correspondingly large uncertainty in the flux estimates 


associated’ with the mean vertical velocity (~ 25%) but they 


still remain significant. 
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TABLE 12 


Sensible Heat Flux Associated with Mean Vertical Wind 


30-Min Means October 21 


30 Mins Ww Heat Flux , Turbulent Heat Total Ho 
Ending (cm/sec) Due w (mW/cm”) Flux (mW/cm*) Flux (mW/cm~*) 
2000 0 4.8 : Sarge dp sa 
2015 5.6 4.0 3 HA lire: 
2030 6.0 4.9 Si26 3.5 
PRIS 9.6 ew) isan je Wee: 
2130 8.8 63 13.0 L923 
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TABLE 13 
Latent Heat Flux Associated With Mean Vertical Wind 


30-Minute Means October 21 


30 Mins W Latent Heat Flux Turbulent Total 
Ending (cm/sec) Due w (mW/cm¢) Latent Heat (mW/cm?) 
Flux (mW/cm2) 


2000 5.8 dees) 4.8 6.6 
2015 ie 0 Lt e8) Dieel Tate 
2030 6.8 18 6.6 8.4 
2125 9.6 PANS. 5.5 Oe 
2130 &.8 2.4 7.6 10.0 
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CHAPTER V 
CONCLUSIONS 


The data collection program at Lake Wabamun in 
1972 was in the nature of a preliminary investigation 
and interpretation of the results must take account of 
the severe limitations of the program in both space (one 
point) and time (parts of three days). Any attempt to 
extend the observed August ratio of fluxes over plume 
and lake to annual, seasonal or even monthly averages 
would be unjustified. Estimates of fluxes from the heated 
water plume have the additional complication of the 
changing size and shape of the plume itself and the 
changing volume of water discharged by the generating 
stations. Conclusions about the enhancement of evaporation 
and heat flux from the lake because of the hot water added 
to it must at present be almost totally qualitative. 

In August, the time of the comparative observations 
of fluxes from the plume and the lake, the difference 
between plume and lake surface temperature should be 
approaching a minimum. Because the fluxes of heat and 
water vapour are directly dependent upon surface temperature, 
the difference in sensible and latent heat fluxes from the 


two surfaces should also be a minimum. Thus the value of 
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approximately 5 (see Table 7) for the ratios of plume 
evaporation and heat flux to lake Lisbbueraton and heat 
flux should be a conservative one with the maximum value 
for this ratio occurring in the winter when all but the 
heated water is frozen. 

If the plumes on the average cover 5% of the lake 
surface (see Nursall et al. 1972) and if the August data 
were representative then it could be concluded that the 
evaporation is enhanced 25% by the heated water. However, 
the question of representative data must remain unanswered 
until INMcre’ data are availables* In sddiicnon this estimate 
in no way takes account of the heat transferred by mixing 
and conduction into the rest of the lake possibly resulting 
ina net rise of temperature in the lake as a whole. The 
additional heat flux may not be significant when considered 
as part of the total energy budget of the lake but the 
increased evaporation may be a more serious problem (see 
Nursall tet) athero 72). 

Because profile measurements were available only 
for October 21, there were not enough data for a significant 
comparison of profile and fast-response estimates of the 
fluxes in general. However, a few comments about the 
comparison for free convection conditions should be valid 


keeping in mind the variation in surface temperature not 


allowed for in profile theory. The two estimates of 
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heat flux were fairly well correlated being generally 
within 15% of one another for Agee ee e means. The estimates 
of water vapour flux for the same intervals were not nearly 
as consistent with each other and in fact the profile 
estimates were greater by approximately 30%. This discrepancy 
may have been related to the near saturation conditions 
and to instrument error because of the presence of liquid 
water in the air but a full explanation is not available. 

The profile estimates of momentum flux were in 
reasonable agreement with the momentum transported by the 
small scale eddies but of course, Dolio dacs was made 
in the profile method for the momentum carried in the large- 
scale convective eddies. The profile procedure may provide 
good estimates of momentum flux in near-neutral stability 
or strong Wind situations but in light wind and free 
convection conditions the picture from gradient information 
is over-simplified. 

The dense fog at the lake in October was found to 
contribute to a 10% increase in the latent heat flux estimates 
but more accurate measurements of this quantity are necessary. 
With lower temperatures the liguid water flux should become 
more important in a fog situation and might have to be 
considered more rigorously. 

The sensible and latent heat fluxes because of the 


non-zero vertical wind were found to contribute significantly 
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to the total heat flux. This indicates that, in the winter 
months when mean upward vertical Pedinie over the heated 
water plumes are likely, the heat flux and evaporation are 
enhanced to a greater degree by the plumes than would be 
expected from the August comparison measurements. 

Finally, the discovery of convective eddies above 
the plume transporting momentum upward against the horizontal 
velocity gradient during light winds was one of the most 
interesting products of the October field trip. These eddies 
were smaller than the heated water in extent but were of 
a larger scale than the ordinary boundary layer turbulence. 
Their occurrence at the tower site appeared to be random 
in nature with an observed duration of not less than 30 
seconds. It was conjectured that these eddies were fairly 
extensive in the vertical with the vertical axis being 
tilted in the downwind direction in such a way that motion 
parallel to this axis would transport momentum upward. No 
information was available on the horizontal extent of the 
eddies and it is not known whether they are symmetrical in 
the cross-wind direction or extend the width of the plume. 
The latter might produce a structure similar to that of the 
large-scale convection forming cumulus in cloud streets 
higher in the troposphere as is frequently observed. With 
the source of these eddies fixed to the heated water plume 


it-should be a relatively easy and rewarding task to gather 
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more information on this convective turbulence at Lake 
Wabamun. 

The problems encountered in this project emphasize 
the necessity of having profile observations outside the 
area of influence of the heated water, especially when large 
temperature differences between the air and water exist. 

As well, the largest part of the temperature gradient was 
found to be below the lowest level of profile measurement. 
More detailed profile observations in the low levels should 
therefore be worthwhile in any future data collection at 


Lake Wabamun. 
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